Qubit residence time measurements with a Bose-Einstein condensate 
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We show that an electrostatic qubit located near a Bose-Einstein condensate trapped in a sym- 
metric double-well potential can be used to measure the duration the qubit has spent in one of its 
quantum states. The strong, medium and weak measurement regimes are analysed and a new type 
of Zeno effect is discussed. The analogy between the residence and the traversal (tunnelling) times 
is highlighted. 
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With the recent progress in quantum information tech- 
nology there often arises a necessity to measure and con- 
trol the state of a two-level quantum system (qubit) . This 
can be achieved by constructing hybrid devices in which a 
microscopic irreversible current between two reservoirs is 
effectively controlled by the qubit's quantum state. Such 
a device can be realised, for example, by placingan elec- 
trostatic qubit close to a point contact (PC) or 
an non-interacting Bose-Einstein condensate trapped in 
a symmetric optical dipole trap '3]. The two systems 
have been shown to affect the observed qubit differently: 
whereas a PC converts any qubit's initial state |i) into a 
statistical mixture exponentially in time Q , decoherence 
of a qubit coupled to a BEC is much slower (~ ^/Vt) 
and strongly dependent on the choice of \i) [3|]. While in 
a PC set up one measures the current across the contact, 
a BEC device is best suited for observing the number 
of atoms which have tunnelled into a previously empty 
reservoir after a time T. The purpose of this Letter is to 
demonstrate that a symmetric BEC device whose Rabi 
oscillations are effectively blocked by the presence of the 
electron in the first qubit's dot, performs a quantum mea- 
surement of the qubit's residence time, i.e., the net du- 
ration the electron has spent in the second dot between 
t = and t = T. We will show that conceptually the 
question of residence time is closely related to the traver- 
sal (tunnelling) time problem still actively debated in the 
literature (see, for example, @]-@) In both cases the time 
in question is the duration a system spends in a specified 
sub-space of its Hilbert space, the sub-barrier region or 
the state in one of the quantum dots. Both quantities 
relate to to the total duration of the system's motion, 
rather than to a single instant, and are conveniently rep- 
resented by a Feynman functional. One can extend von 
Neumann's measurement theory to such functionals [7j, 
but, as far as we know, the BEC device proposed in this 
Letter offers the first practical realisation of such a mea- 
surement. 

The role of a BEC as a measurement tool is best illus- 
trated by considering first somewhat simpler case of a 
condensate coupled to a two-level fluctuator, i.e. a clas- 
sical bistable system switching randomly between two po- 
sitions so that its path q(t) is a random function taking 



values of either or 1 [8|]. Assuming that the tunnelling 
rate of the BEC atoms is enhanced (the barrier is low- 
ered) when g = f , we write the Hamiltonian as 



H 



BEC 



(t) = [n + q{t)sn)]{clcR + c+CL).. sn>o (i) 



The condensate consists of N atoms initially (at t = 0) 
located in the left well. After a time T we wish to count 
the number of atoms in the right well, n, in order to 
obtain information about the noise q{t) for < t < T. 
The probability amplitude for n atoms to tunnel into the 
right well is a functional on the fluctuator's path q given 
by 

Gn^oki-)] = {n\exp[~i{nT + SnT)ic+CR + c+CL)m{2) 

where |n) denotes the BEC state with n atoms in the 
right well, and r is the duration the fluctuator has spent 
in the state q = I, explicitly given by an expression sim- 
ilar to the traversal time functional Tab[x{.)] of Ref. Q 
{6ij is the Kroneker delta) 



Ti[g(.)] = / (5g(t)idt. 



(3) 



For simplicity we will assume that no tunnelling occurs 
for q = 0, i.e. = 0, and the BEC consists of a large 
number of identical non-interacting atoms whose Rabi 
period 2tt/60, is large compared to the observation time 
T, 



N ^oo, sn = a/N^/"^ -> 0, snr o. 



(4) 



Condition (|4]) ensures that if the barrier is permanently 
lowered, {q{t) = 1), there is an irreversible macroscopic 
current into the right reservoir, with the number of tun- 
nelled atoms increasing as a^t^. The spectrum of the 
operator in the exponent of Eq.([5]) consists of equidis- 
tant levels, e„ = (2n — N)5^, n = 0, 1..7V. Expanding 
the exponential in Eq.([2]) in the basis of the correspond- 
ing eigenstates and using the Sterling formula for the 
factorials yields 



Gn^oir) sa Q;"r"exp 



*n>>l 



(-aVV2)/(n!)V2 (5) 
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where 



^ ,^1/2 



la 



(6) 



is the time after which on average n atoms escape to the 
right weU. If the fluctuator's paths are distributed with 
a functional density the probability to find n 

atoms in the right well, |G'„^o(''')P) must be averaged 
further, and we obtain 



P„^o(T)=/ |G„^o(r)pW^(r,r)dr (7) 
Jo 



where the restricted path sum 



W{t,T)^ J2 Sir - rM-)])Wm] 



(8) 



paths 



is the fluctuator residence time probability distribution. 
Thus, finding at t = T exactly n tunnelled atoms allows 
us to conclude that the fluctuator has kept the barrier 
open for a duration r„ — l/\/2a < r < r„ + l/\/2a, i.e. 
that we have measured its residence time to an accuracy 
At = 1/V2a. 

Next we replace the fluctuator with a qubit placed near 
the BEC dipole trap in such a way that the BEC tun- 
nelling rate is enhanced whenever qubit's electron is lo- 
cated in the state The Hamiltonian of the system 
can be written as H = Hq + Hbec, where 



Hq{e) = ea^ai + uo[aia2 + a^ai) 
Hbec ^ {^ + ai5VL){c[cR + CrCl) 



(9) 



5n>o 



and a^2•) the creation operators for the qubit's elec- 
tron electron in the first (tunnelling enhanced) and the 
second (tunnelling suppressed) quantum dot, respec- 
tively. In the following we will put the qubit's Rabi fre- 
quency to unity, w = 1, and re-scale other time and en- 
ergy parameters accordingly. Like a two-state fiuctuator, 
a qubit can alternate between the two states, |1) and |2), 
with the important difference that its trajectory q(t) tak- 
ing the values 1 or 2 is a virtual (Feynman) path. To such 
a path one can assign a probability amplitude <I>[g(t)] 
but not, as above, a probability weight P^[g(t)]. We 
must, therefore, evaluate the number of tunnelled atoms 
at i = T without being able to predict, even with a proba- 
bility, whether the barrier was up or down at any previous 
time < t < T j3] . For a qubiit starting its motion (pre- 
selected) in the state \i) and then at i = T observed (post- 
selected) in a final state |/), this probability amplitude 
is given by ^f^'[q] = (/|9(r))(-ic^)^' exp(-zeT)(g(0)|z), 
where j is the number of times the path crosses from 
one state to another. Following Feynman and Vernon 
[lot we can obtain the probability amplitude ^{^o(T) 
for finding n atoms in the right well given the initial and 
final states of the qubit by multiplying the amplitude in 
Eq.® by $^*~*[(7] and summing over all qubit's paths. 



Assuming, as above, = and recalling that Gnt~o[q{t)] 
only depends on the path's residence time ([3]), we write 



where the restricted path sum (c.f. Refs. 0) 

paths 



(10) 



(11) 



is the qubit's residence time probability amplitude dis- 
tribution. Thus, the quantum analogue of Eq.Q is 



Jo 



(12) 



From Eq.® it is readily seen that the probability pIUq 
results from the interference between the paths with 
T„ — 1/a < r < T„ + 1/a, so that by determining n 
we perform a measurement of the qubit's residence time 
11 1 to a quantum accuracy A'^t = 1/a [13]. Finally, if 
the maximum number of atoms which can tunnel over the 
time T is large, N^ax ~ a^T^ >> 1, we can introduce 
probability density w-^^^ (r, T) for the measured values 



of 



Explicitly we have 



(r,r) = Pl^'oiT)idT,Jdn)-^ = 2a^P, 



n^O- 



W 



(r,r)«(2A)i/2a 



(13) 



I / cxp[-a2(r-T')^]$^^'(T',T)dT'| 
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The measurement statistics are determined by the dis- 
tribution dill) , some of whose properties have been dis- 
cussed in [13|. In particlular, it follows from Eq . lfTTj) that 



(27r)-iexp(-ieT) / exp(ar)(/|C/(r, A)|i)dA 



(14) 



where U{T, A) is the evolution operator for an asymmet- 
ric qubit with the Hamiltonial Hq{X) = Aa^ai -I- (af 02 -I- 
a^fli), whose matrix elements, U^k' = {k\U {T ^ X)\k') are 
given by 

Uii = [cos(£:r/2) - iX£-'^ sin(f T/2)] cxp(-ar/2) 

= exp(— iAT) + uii(A) 
U22 = [cos(£r/2) + iX£-'^ sin(£T/2)] exp(-ar/2) 

= 1+W22(A) (15) 

U12 = -2i£-^sm{£T/2)exp{-i\T/2) ^ U21 

where 5(A) = {X'^+A)^/'^ and 1*11(22) (A) for |A| ^ 00. 
Inserting Eqs. p5)) into Eq. p4|) shows that 



$1^1 (t, T) = 6{t ~T) + (r, T) = $2^2 _ 

$2^1(t,T) = <j)^^\T,T) = ^^^^{t,T) (16) 
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where S{z) is the Dirac delta- function and (j)^^"^ are 
smooth functions of r. For T >> 1 ^-''^^ can be evalu- 
ated by the stationary phase method Considering 
for simplicity a symmetric qubit, e = 0, and introducing 
a new variable ^ = r/T — 1/2, < ^ < 1, wc obtain the 
large-time semiclassical asymptotes valid for < t < T, 

01*-1(t,T) « (2/7rT)i/2 X (17) 
(1 + 2e)i/'*(l - 2^)-"^/^ cos[(l - 4^2)1/2^ + 7r/4] 



-i(2/7rr)i/2(l 



4^^)1/2 gjj^[(^ _ 4^2^i/2y ^ 



(18) 



The oscillatory distributions <^^^^{t,T) and ^^^^^^^y) 
are shown in Fig.l. It is readily seen that after many 
Rabi periods of the qubit, T >> 1, develops a 

stationary region of the width ^ T^/"^ centred at t — 
T/2, which suggests that on average the qubit shares 
its time equally between the states |1) and |2). At the 
same time the singular term 5{t — T) appears to imply 
that the qubit has never left the state |1). There is. 




uT ma 



FIG. 1: Residence time amplitude distributions $^^^(t) and 
$^'~^(r) vs. T and r for 0.02 < r/T < 0.98. Arrows indicate 
the stationary phase region. 

however, no contradiction and next we will show that 
the two conflicting scenarios correspond to two different 
accuracies of the BEC meter and, therefore, are never 
observed at the same time. 

Indeed, for a medium accuracy, 
main contribution to integral (|13p comes from the sta- 



tionary region in Fig.l. and we have 

w'^^{t,T)^w^^^{t,T)= (19) 
(2/^)i/2acos2(T)exph2a2(r-T/2)2] 

(2/^)i/2asin2(r) exp[-2a2(T - T/2f]. 



The Gaussian distributions p9l) . shown in Fig. 2a for 
LuT = 100 and A'r/T = 0.1 by dashed lines, are con- 
sistent with the qubit spending in the state |1) roughly 
half of the total time T. Note that here the contribution 
from the S{t — T) term is cancelled by the oscillations of 
the regular part of <I>^^^(r, T) near t k,T. To model an 
actual measurement and check the accuracy of Eqs. fTH]) 
and (Uni) we have divided the time interval [0, T] into 
Nbin — 100 equal subintervals 5t ~ T /Nun , summed 
the probabilities P^T^l in Ea. p2|) within each interval 
and divided the sum by 5t. The results of this binning 
procedure are shown in Fig. la by the solid lines. As we 




FIG. 2: a) Distributions ■w^'^\t,T) for ojT = 100 and 
At'^ /T = 0.1: Ea. (ll9|l (dashed) and the binning procedure 
with Nhin = 100 (solid); b and c) Same as (a) but for 
ArVr = 0.005. 

increase the coupling strength a, the integral will 
still vanish wherever oscillations of $ are fast compared 
to A'^T = \/a. Where A'^r is small compared to the 
oscillation's period we obtain 



w 



(r,r) w (27r)i/2a-i|$/^'(T,T)|2, < r < r(20) 



Thus, as the accuracy improves, the measurement will 
resolve the pattern of |$^^*(T, t)^ in ever greater de- 
tail. We also note that the probability densities in 



Ea. ((20l) decrease as 



as interaction with BEC sup- 



presses qubit's transitions between the states |1) and |2). 
The approximation (j20p and the results of a binning pro- 
cedure with Nhin ~ 100 are shown in Figs. 2b and 2c for 
IS.'^T jT = 0.005 by the dashed and the solid lines, respec- 
tively. 
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In the high accuracy hmit a ^ the probabihty is con- 
served owing to the S(t — T) and 5[t) terms present, 
as seen from Eqs.dH]), in $i^i(r,T) and <^'^*-'^(T,t), 
respectively. Inserting them into Eq. p^ shows that 
while ^1^2 aj^j ^2^1 (r, T) vanish, w^^i (t, T) and 



'{t,T) become 



linia- 



(r,T) 



(21) 



(2/7r 



,1/2, 



lirrir 



2^2 



;exp[-2a2(r - T)^] w ,5(r - T) 
(T,T)«5(r) 



indicating that the qubit is trapped in its initial state. 
Note that no atoms will tunnel for a qubit starting in 
the second state, |i) = |2), whereas for |i) = |1) one 
obtains a narrow Poisson distribution lima—,ooPn^^ = 
(aT)2"exp(-a2T2)/n!. This variant of Zeno effect, 
which arises not from frequent observations of the mea- 
sured system [l^, but from one accurate evaluation the 
functional ^ over a finite period of time, should be com- 
mon to measurements of other quantities which are non- 
local in time. 

In order to avoid the back action of the BEC on the 
qubit 's evolution and find the 'unperturbed' residence 
time one may be tempted to decrease the coupling by 
putting a ^ 0. Again, it is instructive to analyse 
first the case of a classical fluctuator. In this weak 
couphng hmit Eq.® yields Gn^o{T) ~ a"r"/v^ and 



from Eq.dZD we obtain Pn^^iT) 



^,2n/^2n 



) /nl where 



2n\ 



/p T^'^W{T)dT is the n-th even moment of the 



probability distribution W{t) > 0. Thus from the ratio 
Pi^o(r)/Po^o(r) w a'^ir^) we can determine (t^) and, 
should the dispersion be small, the mean residence time 
(r) « yp). 

Similarly, for a qubit from Eqs.©, and in 

the limit a -> we find Pl^' / pI'^' « a2n|^-„|2/^l 
with ^ = /r"$^"'*(T,T)dT//$/-*(T,r)dr = 
{—i)^d\log{J\U{T,X)\i)\\=Q. In particular, we have 



Pit 



-0 / ^0+- 



2 I— 12 



\fY ^a^\T/2 + t'An{T)/2Y. (22) 



where r is the weak value of the residence time analogous 
to the Larmor tunnellin g ti me first introduced to quan- 
tum scattering by Baz' [l6|. It diverges whenever Rabi 
oscillations put the unperturbed qubit into the state |2), 
T = {k + l/2)7r, fc = 0, 1..., may exceed the total dura- 
tion of motion T, and cannot be interpreted as a valid 
residence time. This problem is common to all weak mea- 
surements introduced in [l7|, whose accuracy is so poor 
that they do not destroy coherence between different val- 
ues of the measured quantity [l^. The weak residence 
time f in Eq. ((22|) is the first moment of an alternating 
amplitude distribution ^^^^(T,t) and as such is not di- 
rectly linked to the physical values < t < T [ist - 
In summary, we have shown that a hybrid device consist- 
ing of an electrostatic qubit coupled to a BEC trapped 
in a symmetric double-well potential can be used to per- 



form the qubit's residence time measurements. Depend- 
ing on the strength of the coupling the measurement can 
be 'weak' or strong. An accurate (strong) measurement 
leads to a finite time Zeno effect trapping the qubit's elec- 
tron in one of the quantum dots. Mathematical explana- 
tion linking the effect to the presence of singular terms 
in the residence time amplitude distribution should also 
apply to a wide range of similar measurements. 
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